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Abstract 

The constraints imposed on the Minimal Supersymmetric Standard Model (MSSM) parameter 
space by the Large Hadron Collider (LHC) Higgs mass limit and gluino mass lower bound are 
revisited. We also analyze the thermal relic abundance of lightest neutralino, which is the Lightest 
Supersymmetric Particle (LSP). We show that the combined LHC and relic abundance constraints 
rule out most of the MSSM parameter space except a very narrow region with very large tan /? (~ 
50). Within this region, we emphasize that the spin-independent scattering cross section of the 
LSP with a proton is less than the latest Large Underground Xenon (LUX) limit by at least 
two order of magnitudes. Finally, we argue that non-thermal Dark Matter (DM) scenario may 
relax the constraints imposed on the MSSM parameter space. Namely, the following regions are 
obtained; mo — 0(4) TeV and mi /2 — 600 GeV for low tan/3 (~ 10); mo ~ rni /2 — 0(1) TeV or 
mo ^ 0(4) TeV and mi /2 — 700 GeV for large tan/3 (~ 50). 
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I. INTRODUCTION 


The most recent observations by the Planck satellite confirmed that 26.8% of the universe 
content in the form of DM and the usual visible matter only accounts for 5% [T]. The LSP 
remains one of the best candidates for the DM [a [3]. It is a Weakly Interacting Massive 
Particle (WIMP) that can naturally account for the observed relic density of DM. 

Despite the absence of direct experimental verihcation, Supersymmetry (SUSY) is still 
the most promising candidate for a unified theory beyond the Standard Model (SM). SUSY 
is a generalization of the space-time symmetries of the quantum held theory that links the 
matter particles (quarks and leptons) with the force-carrying particles, and implies that 
there are additional ‘superparticles’ necessary to complete the symmetry. In this regards, 
SUSY solves the problem of the quadratic divergence in the Higgs sector of the SM in a very 
elegant natural way. The most simple supersymmetric extension of the SM, which is the 
most widely studied, is know as the MSSM nn. In this model, certain universality of soft 
SUSY breaking terms is assumed at grand unihcation scale. Therefore, the SUSY spectrum 
is determined by the following four parameters: universal scalar mass mo, universal gaugino 
mass mi/ 2 , universal trilinear coupling Aq, the ratio of the vacuum expectation values of 
Higgs bosons tan f3. In addition, due to i?-parity conservation, SUSY particles are produced 
or destroyed only in pairs and therefore the LSP is absolutely stable, implying that it might 
constitute a possible candidate for DM, as first suggested by Goldberg in 1983 [7]. So 
although the original motivation of SUSY has nothing to do with the DM problem, it turns 
out that it provides a stable neutral particle and, hence, a candidate for solving the DM 
problem. 

The landmark discovery of the SM-like Higgs boson at the LHC, with mass ~ 125 GeV 
[8], might be an indication for the presence of SUSY. Indeed, the MSSM predicts that there 
is an upper bound of 130 GeV on the Higgs mass. However, this mass of lightest Higgs boson 
implies that the SUSY particles are quite heavy. This may justify the negative searches for 
SUSY at the LHC-run I [9]. However, it is clearly generating a new ‘little hierarchy problem’. 

Moreover, the relic density data [1] and upper limits on the DM scattering cross sections 
on nuclei (LUX 110] and other direct detection experiments [H]) impose stringent constraints 
on the parameter space of the MSSM |T2]. In fact, combining the collider, astrophysics and 
rare decay constraints [T3HTH] almost rule out the MSSM. It is tempting therefore to explore 
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well motivated extensions of the MSSM, such as NMSSM [19] and BLSSM [20], which may 
alleviate the little hierarchy problem of the MSSM through additional contributions to Higgs 
mass [m 121] and also provide new DM candidates [22] that may account for the relic density 
with no conflict with other phenomenological constraints. 

In this article we analyze the constraints imposed by the Higgs mass limit and the gluino 
lower bound, which are the most stringent collider constraints, on the MSSM parameter 
space. In particular, these constraints imply that the gaugino mass, mi/ 2 , resides within 
the mass range: 620 GeV < mi /2 ^ 2000 GeV. While the other parameters are much less 
constrained. We study the effect of the measured DM relic density on the MSSM allowed 
parameter space. We emphasized that in this case all parameter space is ruled out except 
few points around tan/? ~ 50, mo ~ 1 TeV and mi /2 ~ 1.5 TeV. We also investigate the 
direct detection rate of the LSP at these allowed points in light of the latest LUX result. 
Finally we show that if one assumes non-standard scenario of cosmology with low reheating 
temperature, where the LSP may reach equilibrium before the reheating time, then the relic 
abundance constraints on (mo,mi/ 2 ) can be signihcantly relaxed. 

The paper is organized as follows. In section 2 we briefly introduce the MSSM and study 
the constraints on the (mo,mi/ 2 ) plane from Higgs and gluino mass experimental limits. In 
section 3 we study the thermal relic abundance of the LSP in the allowed region of parameter 
space. We show that the combined LHG and relic abundance constraints rule out most of 
the parameter space except the case of very large tan (3. We also provide the expected rate of 
direct LSP detection at these points with large tan [3 and TeV masses. Section 4 is devoted 
for non-thermal scenario of DM and how it can relax the constraints imposed on MSSM 
parameter space. Finally we give our conclusions in section 5. 


II. MSSM AFTER THE LHC RUN-I 

The particle content of the MSSM is three generations of (chiral) quark and lepton super- 
helds, the (vector) superhelds necessary to gauge the SU (3)c x SU{2)l x U (l)y gauge of the 
SM, and two (chiral) SU{2) doublet Higgs superhelds. The introduction of a second Higgs 
doublet is necessary in order to cancel the anomalies produced by the fermionic members of 
the hrst Higgs superheld, and also to give masses to both up and down type quarks. The 
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interactions between Higgs and matter superfields are described by the superpotential 

W = huQLUlH2 + hDQlDlHi + hLlLElHi + /ih/ih/2. (1) 

Here Ql contains SU{2) (s)quark doublets and t/£, D^i are the corresponding singlets, 
(s)lepton doublets and singlets reside in and i?£ respectively. While Hi and H 2 denote 
Higgs superfields with hypercharge H = =f|. Further, due to the fact that Higgs and lepton 
doublet superfields have the same SU{3) x SU{2)l x U{1)y quantum numbers, we have 
additional terms that can be written as 


W' = X^jkL.L.El + KjkLiQiDl + KjkDtD^^ + /i.L.h/s- 


( 2 ) 


These terms violate baryon and lepton number explicitly and lead to proton decay at unac¬ 
ceptable rates. To forbid these terms a new symmetry, called i?-parity, is introduced, which 
is defined as Rp = (^—where B and L are baryon and lepton number and S is the 
spin. There are two remarkable phenomenological implications of the presence of i?-parity: 
i) SUSY particles are produced or destroyed only in pair, ii) The LSP is absolutely stable 
and, hence, it might constitute a possible candidate for DM. 

In MSSM, a certain universality of soft SUSY breaking terms at grand unification scale 
Mx = 3 X 10^® GeV is assumed. These terms are defined as mo, the universal scalar soft 
mass, mi/ 2 , the universal gaugino mass, Aq, the universal trilinear coupling, B, the bilinear 
coupling (the soft mixing between the Higgs scalars). In order to discuss the physical 
implication of soft SUSY breaking at low energy, we need to renormalize these parameters 
from Mx down to electroweak scale. In addition the MSSM contains another two free SUSY 
parameters: /i and tan/3 = {H 2 )/{Hi). Two of these free parameters, /i and B, can be 
determined by the electroweak breaking conditions: 

^H2 


2 

h = 


mp-„ tan^ (3 


tan^ /3 — 1 


- Mi/2, 


. -2m| 

sm2/3 = 


m\ + m 2 


(3) 

(4) 


Thus, the MSSM has only four independent free parameters: mg, mi/ 2 , Hq, tan/3, besides to 
the sign of /r, that determine the whole spectrum. 

In the MSSM, the mass of the lightest Higgs state can be approximated, at the one-loop 
level, as [23] 


mf^<Mz + 


m: 


sin^ /3 


log 


9 2 

m~ m~ 


ml 


(5) 
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FIG. 1: MSSM parameter space for tan/3 = 10 (top panel) and 50 (bottom panel) with Aq = 0 
and 2 TeV. The green region indicates for 124 < rrih < 126 GeV. The blue region is excluded 
because the lightest neutralino is not the LSP. The pink region is excluded due to absence of 
radiative electroweak symmetry breaking {fj? becomes negative). The gray shadow lines denote 
the excluded area because of nig <1.4 TeV. 

Therefore, if one assumes that the stop masses are of order TeV, then the one-loop effect 
leads to a correction of order 0(100) GeV, which implies that 

< v^(90 GeV)2 + (100 GeV)^ ~ 135 GeV. (6) 

The two-loop corrections reduce this upper bound by a few GeVs |2^. Hence, the MSSM 
predicts the following upper bound for the Higgs mass: nih < 130 GeV, which was consistent 
with the measured value of Higgs mass (of order 125 GeV) at the LHG |8]. 

In Fig. we display the contour plot of the SM-like Higgs boson: nih € [124,126] GeV 
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in (mo,mi/ 2 ) plane for different values of Aq and tan/3. It is remarkable that the smaller 
Aq is, the smaller mi /2 is needed to satisfy this value of Higgs mass. It is also clear that 
the scalar mass itiq remains essentially unconstrained by Higgs mass limit. It can vary from 
few hundred GeVs to few TeVs. Such large values of mi /2 seem to imply a quite heavy 
SUSY spectrum, much heavier that the lower bound imposed by direct searches at the LHC 
experiments in centre of mass energies i/s = 7, 8 TeV and total integrated luminosity of 
order 20 fb“^. Furthermore, the LHC lower limit on the gluino mass: rrig > 1.4 TeV |25] . 
excluded the values of mi /2 < 620 GeV that was allowed by Higgs mass constraints for 
mo > 4 TeV. Furthermore, this region is shown with dashed lines in Fig. 


III. DARK MATTER CONSTRAINTS ON MSSM PARAMETER SPACE 


A. The LSP as dark matter candidate 


The neutralinos Xi (1=1)2,3,4) are the physical (mass) superpositions of two fermionic 
partners of the two neutral gauge bosons, called gaugino (bino) and (wino), and of 
the two neutral Higgs bosons, called Higgsinos and H^- The neutralino mass matrix is 
given by [2B] 


/ 


Ml 


0 


—Mz cos (3 sin 6w 


Mz sin (3 sin 6w 


\ 


Mn = 


0 


M2 


Mz cos f3 cos 0w —Mz sin (3 cos 6w 


(7) 


—Mz cos j3 sin 9w Mz cos j3 cos 9\y 0 —p 

^ Mzsml3sm9w — Mz sin(3 cos 9w —P 0 j 

where Mi and M 2 are related due to the universality of the gaugino masses at the grand 
unihcation scale. Mi = M 2 , where pi, p 2 are the gauge couplings of U{1)y and SU{2)l 

^92 

respectively. This Hermitian matrix is diagonlaized by a unitary transformation of the 
neutralino helds, M^“^ = N'^M^N. The lightest eigenvalue of this matrix and the corre¬ 
sponding eigenstate say y has good chance of being the LSP. The lightest neutralino will be 
a linear combination of the original helds: 


X = NuB° + Yi2lU° + + NuHl 


( 8 ) 


The phenomenology and cosmology of the neutralino are governed primarily by its mass and 
composition. A useful parameter for describing the neutralino composition is the gaugino 
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FIG. 2: The mass of lightest neutralino versus the purity function in the region of parameter space 
allowed by gluino and Higgs mass limits. 

“purity” function fg = \Nn\^ + \Nu\'^ |26]. If/g > 0.5, then the neutralino is primarily gaug- 
ino and if fg < 0.5, then the neutralino is primarily Higgsino. Actually if |/i| > IM 2 I > Mz, 
the two lightest neutralino states will be determined by the gaugino components, similarly, 
the light chargino will be mostly a charged wino. While if |/r| < IM 2 I, the two lighter neu- 
tralinos and the lighter chargino are all mostly Higgsinos, with mass close to |p|. Finally if 
|/i| ~ \M 2 \, the states will be strongly mixed. 

Here, two remarks are in order: i) The above mentioned constraints in mi /2 from Higgs 
mass limit and gluino mass lower bound imply that > 240 GeV, which is larger than the 
limits obtained from direct searches at the LHC. Moreover, an upper bound of order one 
TeV is also obtained (from Higgs mass constraint), ii) In this region of allowed parameter 
space, the LSP is essentially pure bino, as shown in Fig. This can be easily understood 
from the fact that /i-parameter, determined by the radiative electroweak breaking condition, 
Eq. ([^, is typically of of order mo and hence it is much heavies than the gaugino mass Mi. 

B. Relic denisty 

As advocated in the previous section, the LSP in MSSM, the lightest neutralino y, is 
a perfect candidate for DM. Here, we assume that y was in thermal equilibrium with the 
SM particles in the early universe and decoupled when it was non-relativistic. Once the 
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X annihilation rate '^)^x dropped below the expansion rate of the universe, 

< H, the LSP particles stop to annihilate, fall out of equilibrium and their relic density 
remains intact till now. The above v) refers to thermally averaged total cross section 

for annihilation of yy into lighter particles times the relative velocity, v. 

The relic density is then determined by the Boltzmann equation for the LSP number 
density (n^) and the law of entropy conservation: 


dn^ 

dt 

ds 

dt 


-3Hn^ - v) [(n,)^ - (n<^f] , 
—3Hs, 


(9) 

( 10 ) 


where n'^ is the LSP equilibrium number density which, as function of temperature T, 
is given by Here and are the mass and the number 

of degrees of freedom of the LSP respectively. Finally, s is the entropy density. In the 
standard cosmology, the Hubble parameter H is given by H{T) = where Mpi = 

1.22 X 10^® GeV and is the number of relativistic degrees of freedom. Let us introduce 
the variable x = rn^/T and define Y = n^/s with Yeq = s. In this case, the Boltzmann 
equation is given by 


dY 1 ds 
dx 3Hdx^^ 

In radiation domination era, the entropy, as function of the temperature, is given by 

27r2 

s(i) = X 


( 11 ) 


( 12 ) 


which is deduced from the fact that s = (p + p)/T and g^^ is the effective degrees of freedom 
for the entropy density. Therefore one finds 


ds 3s 

dx X 


( 13 ) 


Thus, the following expression for the Boltzmann equation for the LSP number density is 
obtained 


dY (<“” 


v) 


(y^ - Yel) 


(14) 


If one considers the s-wave and p-wave annihilation processes only, the thermal average 


(^ann then shows as 




\ 

v) =a^ + -, 


X 


( 15 ) 




FIG. 3: Feynman diagrams contributing to early-universe neutralino (y^) annihilation into fermions 
through sfermions, Z-gauge boson, and Higgs. 


where and are the s-wave and p-wave contribntions of annihilation processes, respec¬ 
tively. The relic density of the DM candidate is given by 


^1^2 ^ m^soY^{oo) 


Pc/h‘^ 


(16) 


where Sq = 2282.15 x 10 GeV^, Pc = 8.0992 x 10 GeV^, and by solving the Boltzmann 
eqnation, one can hnd D^(c)o) as follows 1271 

-1 


+ 


36. 


A, \x(Tf) ' xHT,) 


( 17 ) 


where Ty is the freeze-ont temperature, = s{m-^/Hiprriy^) and x{Tf) is given by 


xiTf) = In 


ci^A^c(c -|- 2) 


+ 


_^A 

<Tf)J 


(18) 


where value c = | results in a typical accuracy of about 5 — 10% 

more than sufficient for our purposes here. 

The lightest neutralino may annihilate into fermion-antifermion (//), W^W~, ZZ, 
W^H~, ZA^ ZH^ Zh^ H^H~ and all other contribution of neutral Higgs. For a bino- 
like LSP, i.e. Nn ~ 1 and Nn ~ 0, i = 2,3,4, one finds that the relevant annihilation 
channels are the fermion-antifermion ones, as shown in Fig. and all other channels are 
instead suppressed. Also, the annihilation process mediated by Z gauge boson is suppressed 
due to the small Zxx coupling cx except at the resonance when 

which is no longer possible due to the above mentioned constraints. Furthermore, one finds 
that the annihilation is predominantly into leptons through the exchanges of the three slep- 
ton families (IlJr), with I = e,/i, r. The squarks exchanges are suppressed due to their 
large masses. 
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FIG. 4: LSP relic abundance constraints (red regions) on (ttiq — ^ 21 / 2 ) plane for tan/3 and ^0 a-s 
in Fig. The LUX result is satisfied by the yellow region. The other color codes are as in Fig. 


In Fig.|^we display the constraint from the observed limits of on the plane (mo—mi/ 2 ) 
for Ao = 0, 2000 GeV, tan /3 = 10, 50 and p > 0. Here we used micrOMEGAs [28] to compute 
the complete relic abundance of the lightest neutralino, taking into account the possibility of 
having co-annihilation with the next-to-lightest supersymmetric particle, which is typically 
the lightest stau. In this figure the red regions correspond to a relic abundance within the 
measured limits [1]: 

0.09 < nh^ < 0.14 (19) 

It is noticeable that with low tan (3 (~ 10), this region corresponds to light mi /2 (< 500 GeV), 
where a significant co-annihilation between the LSP and stau took place. However, this 
possibility is now excluded by the Higgs and gluino mass constraints [22] • At large tan/3, 
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FIG. 5: The relic abundance versus the mass of the LSP for different values of tan/3. Red points 
indicate for 40 < tan/3 < 50 and blue points for 30 < tan/3 < 40. All points satisfy the above 
mentioned constraints. 


another region is allowed dne to a possible resonance dne to s-channel annihilation of the 
DM pair into ferniion-antiferniion via the psendoscalar Higgs boson A at Ma — 2m^ EDI- 
For Ho = 0, a very small part of this region is allowed by the Higgs mass constraint, while 
for large Aq (~ 2 TeV) a slight enhancement of this part can be achieved. In Fig. we 
zoom in on this region to show the explicit dependence of the relic abnndance on the LSP 
mass and large valnes of tan /3. As can be seen from this hgnre, there is no point can satisfy 
the relic abnndance stringent constraints with tan /3 < 30. 

C. Direct Detection 

Perhaps the most natnral way of searching for the nentralino DM is provided by direct 
experiments, where the effects indnced in appropriate detectors by nentrali-nonnclens elastic 
scattering may be measnred. The elastic-scattering cross section of the LSP with a given 
nnclens has two contribntions: spin-dependent contribntion arising from Z and q exchange 
diagrams, and spin-independent (scalar) contribntion dne to the Higgs and sqnark exchange 
diagrams, which is typically snppressed. The effective scalar interaction of nentralino with 
a qnark is given by 

^ = aqXXm, ( 20 ) 
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FIG. 6: Spin-independent scattering cross section of the LSP with a proton versus the mass of the 
LSP within the region allowed by all constraints (from the LHC and relic abundance). 

where is the neutralino-quark effective coupling. The scalar cross section of the neutralino 
scattering with target nucleus is given by [2] 

Arry^ 2 

asi=^{Zf, + {A-Z)U)\ ( 21 ) 

71 

where Z and A — Z are the usual atomic numbers, is the reduced mass of the nucleon 
and fp, fn are the neutralino coupling to protons and neutrons respectively. 

In Fig. we display the MSSM prediction for spin-independent scattering cross section 
of the LSP with a proton after imposing the LHC and relic abundance constraints. It is 
clear that our results for are less than the recent LUX bound (blue curve) by at least 
two order of magnitudes. This would explain the negative results of direct searches so far. 

IV. NON-THERMAL DARK MATTER AND MSSM PARAMETER SPACE 

In the previous section, we assumed standard cosmology scenario where the reheating 
temperature Tjih is very large, namely Tjih » Tf ~ 10 GeV. However, the only constraint 
on the reheating temperature, which could be associated with decay of any scalar field, 0, 
not only the inflaton held, is Trh ^ 1 MeV in order not to spoil the successful predictions 
of big bang nucleosynthesis. 

A detailed analysis of the relic density with a low reheating temperature has been 
carried out in Ref. |3Tj. It was emphasized that for a large annihilation cross section. 
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FIG. 7: LSP non-thermal relic abundance constraints (red regions) on (mo — ^ 2 - 2 / 2 ) plane for tan (3 
and ^0 as in Fig. The color codes are as in Fig. 


^^ann > 10“^^ GeV“^ SO that the neutralino reaches equilibrium before reheating, and if 
there is a large number of neutralinos produced by the scalar field 0 decay, then the relic 
density is estimated as [32] 


= 


3 




2 (27r2/45) v) pdS q' 

Here the reheating temperature is dehned as 1271 


Trh = 


90 




1/4 




PI 


, 1/2 


where the decay width F^ is given by 


27r 


( 22 ) 


(23) 


( 24 ) 
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The scale A is the effective suppression scale, which is of order the grand unihcation scale 
Mx- Therefore, for scalar held with mass ~ 10^ GeV one hnds T,^ ~ 10“^^ GeV and 
in our calculations, we have used = 10.75 due to the consideration of a low reheating 
temperature scenario. 

In Fig. we show the constraints imposed on the MSSM (mg — mi/ 2 ) plane in case of 
non-thermal relic abundance of the LSP for tan/3 = 10, 50 and Aq = 0, 2 TeV. In this plot, 
we also imposed the LHG constraints, namely the Higgs mass limit and the gluino mass 
lower bound, similar to the case of thermal scenario. It is clear from this hgure that the 
stringent constraints imposed of the MSSM parameter space by thermal relic abundance are 
now relaxed and now low tan/3 (~ 10) is allowed but with very heavy mo (~ Cl(4) TeV) 
and mi /2 ~ 600 GeV. In addition, the following two regions are now allowed with large 
tan/3 (~ 50): (/) mo ~ mi /2 ~ 0(1) TeV; {ii) mo ^ 0(4) TeV and mi /2 ^ 700 GeV. 
The SUSY spectrum associated with these regions of parameters space could be striking 
signature for non-thermal scenario at the LHG. 


V. CONCLUSION 

We have studied the constraints imposed on the MSSM parameter space by the Higgs 
mass limit and the gluino lower bound, which are the most stringent collider constraints 
obtained from the LHC run-I at energy 8 TeV. We showed that mi /2 resides within the mass 
range: 620 GeV < mi /2 ^ 2000 GeV, while the other parameters (mo, Aq, tan/3) are much 
less constrained. We also studied the effect of the measured DM relic density on the MSSM 
allowed parameter space. It turns out that most of the MSSM parameter space is ruled out 
except few points around tan/3 ~ 50, uiq ~ 1 TeV and mi /2 ~ 1-5 TeV. We calculated the 
spin-independent scattering cross section of the LSP with a proton in this allowed region. 
We showed that our prediction for a'gj is less than the recent LUX bound by at least two 
order of magnitudes. We have also analyzed the non-thermal DM scenario for the LSP. We 
showed that the constraints imposed on the MSSM parameter space is relaxed and low tan /3 
is now allowed with mo — 0(4) TeV and mi /2 — 600 GeV. Also two allowed regions are now 
associated with large tan/3 (~ 50), namely: mg ~ mi /2 — 0(1) TeV or mg — 0(4) TeV and 
mi /2 ~ 700 GeV. 
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